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Abstract 

We present a new cubic theory of gravity in five dimensions which has second order traced field 
equations, analogous to BHT new massive gravity in three dimensions. Moreover, for static spheri- 
cally symmetric spacetimes all the field equations are of second order, and the theory admits a new 
asymptotically locally flat black hole. Furthermore, we prove the uniqueness of this solution, study 
its thermodynamical properties, and show the existence of a C-function for the theory following the 
arguments of Anber and Kastor (arXiv:0802.1290 [hep-th]) in pure Lovelock theories. Finally, we 
include the Einstein-Gauss-Bonnet and cosmological terms and we find new asymptotically AdS 
black holes at the point where the three maximally symmetric solutions of the theory coincide. 
These black holes may also possess a Cauchy horizon. 



Contents 

1 A cubic theory in five dimensions 

2 Field equations for static spherically symmetric spacetimes 

3 Birkhoff's theorem 

4 Entropy function 

5 Black hole thermodynamics 

6 Nonhomogeneous combinations: Asymptotically AdS black holes 

7 Generalization to arbitrary higher order 

8 Further comments 



1 



A General static, spherically symmetric solution for the 
non-homogeneous combination 

B A new quartic theory of gravity in seven dimensions 

1 A cubic theory in five dimensions 

There has been a considerable interest in higher curvature theories of gravity in the last few decades. 
Among them, the most prominent one is Lovelock theory of gravity, which is a natural generalization 
of Einstein's General Relativity in higher dimensions |TJ. Lovelock theories are characterized by the 
special property that the field equations are of at most second order in derivatives of the metric. As 
a consequence, they can have stable, ghost-free, constant curvature vacua [2]. In five dimensions, the 
most general Lovelock gravity is given by an arbitrary linear combination of the cosmological constant, 
the Ricci scalar curvature and the four-dimensional Euler density (also known as the Gauss-Bonnet 
term) which is quadratic in curvature. The Lovelock terms which are higher order in curvature vanish 
identically in five dimensions. However, recently another higher curvature theory in three dimensions 
has drawn a lot of attention. This theory, known as the BHT New Massive Gravity [3J, supplements 
to the usual Einstein-Hilbert term, a precise combination of quadratic curvature invariants. One of 
the key properties of the theory is that the trace of the field equations arising from the pure quadratic 
part, being proportional to itself, is of second order. The pure quadratic part is the unique quadratic 
curvature invariant which possesses this property in three dimensions [3], [5]. In five dimensions, there 
are two linearly independent cubic curvature invariants which share this property [6]. One of them can 
be expressed as a complete contraction of three conformal Weyl tensors. In this section, we present 
another linearly independent cubic curvature invariant (in five dimensions), which shares this property 
with the pure quadratic BHT (in three dimensions). Consider the following action in five dimensions 



d b x , (1) 
where [n] = [Length] is assumed to be positive hereafter, and the Lagrangian is given by 

r T>o.b jjce jjdf p cdn be r>a i ^ Z?a r>b r>c ^DD« D* i ^ E>3 

L ~~Q K cd K bf K ae ~ K ab K cd K e ~ ^ab K c K d + ^ K b K c K a ~ 1^ KK b K a + ■ \ l ) 

Varying the action with respect to the metric gives the following fourth order field equations 

E a b = - - 
6 



3R ahd R b PrdR pgr k ~ ^P^q( RP g q h R ab ~ ^hbg 11 ^) ~ ^abR^cd^nf^ me 



RacbdR CSPq Rpqs ~ lR c db iR qh + R b ^ R adc ' R qli ~ V p V g (R a hR^ hq + R a hR^ P + RbhRa h<1 

(3) 
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+R q h R a \ P + R\K\ h + \{9 pq Ra hcd R bhcd + 9a b R prcd R\ cd - g a p R b rcd R q rcd - 9 b p R a rcd R\ cd )) 



— n £> cdr> ne r>m 
' 2" a b rnn JT cd iL e 



R ac R b fcd R fd + 2R acbd R^ dq R fg + V p V q (R ab R pq - R p R q + g pq R acbd R cd 



+9abR PCqd Rcd ~ 9a P R 9 cbdR Cd ~ 9 b P R 9 ca dR Cd ) ~ iflabRmn R^R^d 



1 

+ 3 
1 
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3R acbd R ec R d + ^V p V,(^i? a c i? fcc + g ab Bf*R t ? - g b p R qc R ac - 9 a p R qc R b c) - \g a bR m n R n c R C r, 
R ab R cd Rcd + 2RR cd R acbd + V p V q (g ab g pq R cd R cd + <f iiiU - g p g b q R cd R c d + g a bRR pq 



-g p RR q - g p RR b q ) - \g a bRR m n R r > 



1 

+ 12 



Si? 2 ^ + 3V p V q (g ab g pq R 2 - g p g b q R 2 ) - ^g ab R 3 



(4) 



Interestingly, the trace of the field equations, being proportional to the Lagrangian ([2]), is of second 
order. Indeed 

E a a = \C . (5) 

The theory defined by action ([T]) has several other interesting aspects which we will discuss in 
the following sections. In section [21 we show that for static spherically symmetric ansatz, the field 
equations reduce to second order. We then integrate the field equations to obtain the most general 
solution in this family, which, for a certain range of the integration constant, describes a (topological) 
black hole. In section [3l we prove a Birkhoff 's theorem for the black hole solution. In section [H 
we show the existence of a C-function for the theory using Wald's formula. We then briefly discuss 
the thermodynamical properties of the topological black hole in section [5j In section [6l we add an 
Einstein-Gauss-Bonnet and cosmological term to the cubic Lagrangian and fix the coupling constants 
so as to have a unique maximally symmetric vacuum. At this special point, we obtain a asymptotically 
AdS black hole solution. Finally, in section [5J we propose a conjecture for the exitance of arbitrary 
higher order Lagrangians, sharing the above features and offer some comments. Appendix A contains 
the general static spherically symmetric solution for the non-homogeneous cubic combinations, and in 
Appendix B we present the quartic generalization in seven dimensions. 



2 Field equations for static spherically symmetric spacetimes 

Let us consider the metric 

ds 2 = -f (r) dt 2 + ^- ) + r 2 dS| , (6) 

where dY>3 is the line element of a Euclidean three-dimensional space of constant curvature 7 = ±1, 0. 
For 7 = +1, since S3 is locally isomorphic to S 3 , the metric © possesses spherical symmetry. In 
the case 7 = — 1, the metric S3 corresponds to an identification of the hyperbolic space H3, while for 
7 = 0, S3 is locally flat. Note that the metric © is the most general static metric which is compatible 



3 



with the isometries of S3. 

The field equations (HJ), when evaluated on the metric ([HD, reduce to 



Et t = ^(9-7) 2 {2g-3g'r-2 1 ) , (7) 

E "r = (5 " 7) 2 (2/5 - 3r<//' - 2 7 /) , (8) 

£i i = 7^1 (7 - <?) f' 2 + ff'r {5g'rg + A 7 g - 4g 2 - jrg') (9) 

+4/ 2 r ( 7 -g)g'- 2fgr 2 (7 - g) f" + 4/ 2 (7 - gf] ^ (10) 

where (') denotes derivative with respect to r, and the indices i,j run along the base manifold £3. 
Then, the field equations for the metric ([B]) reduce to a set of nonlinear second order equations for 
the functions / (r) and g (r). Let us first analyze the nontrivial branch, where g(r) 7^ 7. Solving this 
system, it is easy to see that the only nontrivial solution within the static family (jBJ is 

V / fyf^ji _|_ ry 



r 



where c is an integration constant. This metric is asymptotically locally flat since R a g — > when 
goes to infinity, and has a curvature singularity at the origin, which can be realized by evaluating the 
Ricci scalar 

88c 

*=-^75- < 12 > 

This singularity is hidden by an event horizon when c is positive and 7 = — 1, in which case it is more 
convenient to rewrite the metric (jlip as 

&2 = -((£f-^* 2 + 77fc +r2dEi ' (13> 



r+ 



where r + = c~ 3 / 2 is the location of the horizon, whose geometry is given by H3/T, where T is a freely 
acting, discrete subgroup of O (3, 1). 

In the next section, we prove a Birkhoff 's theorem for the solution when the staticity condition is 
removed. 

Note that there is another branch of solutions for the system ([7|)- (|10p . for which g(r) = 7 and 
f(r) is an arbitrary function. This "degenerated" behavior is also common to all Lovelock theories 
possessing a unique maximally symmetric solution [7]-|12j. 



3 Birkhoff's theorem 

Let us consider the metric 

ds 2 = -f (t, r) dt 2 + -— + r 2 dX 2 , (14) 
9(t,r) 
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where again the manifold S3 is a compact manifold of constant curvaturd^ 7 = ±1,0. For this family 
of spacetimes, the nonvanishing diagonal components of the field equations are 



E t t = ^{9-l?{2g-2,g'r-2 1 ) , (15) 
1 



E r r = 2f^(9~ 1? (2/5 - 3rgf - 2jf) , (16) 

E 'j = JTg^s ^ V ft " 9) f' 2 + ff'9 2 r {5g'rg + 4 7 g - Ag 2 - jrg') (17) 

+4/V r ( 7 -g)g'- 2fg\ 2 (7 - g) f" + 4/V (7 - gf (18) 

+gr 2 (7 ~g)fg + fr 2 (3 7 + g) g 2 - 2fgr 2 (7 - g) g\ 5^ (19) 

where (') and (•) denote partial derivatives with respect to r and t respectively. The off-diagonal 
component E t r contains mixed partial derivatives of the metric functions. Note that equations ()15p 
and (fTUj) have the same expression as their static counterpart ([7]) and (jSJ) respectively, since they do 
not contain derivatives with respect to time. 
Equation (fT5]h is solved by 

g(t,r) =F(t)r 2 / 3 + 7 , (20) 

where F (t) is an arbitrary integration functiorj 2 ]. Inserting this expression for g (t, r) in (I16p we obtain 
the following equation for / (t, r): 

3r (7 + F (t) r 2 / 3 ) - 2/ (t, r) F (t) r 2 / 3 = , (21) 

whose solution is 

f(t,r)=H(t)g(t,r) , (22) 

where H (t) being an integration function, can be absorbed by a time rescaling. Thus, without any 
loss of generality, the equations E\ = and E r r = are solved by 

f(t,r) =g(t,r) = F (t) r 2 / 3 + 7 . (23) 

Now, after inserting (I23P in (I19p . we obtain the following equation for F (i): 

F 2p r 2/3 + ( 2 p2 + 7 = o _ (24) 



This implies that , for 7 7^ 0, F (£) must be a constant c, and the metric (|14p reduces to (jlip . which is 
the static metric obtained previously. Thus, for 7 7^ we have proved the Birkhoff's theorem, since we 
have shown that the most general solution of the theory, within the family of spacetimes f)14[) , is static 
and is given by (jlip . For 7 = 0, equation ()24|) implies that F (t) = et + c (e and c being integration 



In general, in front of £3 there can be a generic function F (t, r) which after a gauge fixing, depending on the norm 
of its gradient can be chosen to be r 2 , t 2 or a constant. 

2 When F(t) is identically vanishing, the remaining field equations are automatically solved for an arbitrary f(t, r) — 
f(r), and we get back the degenerate solution mentioned previously. 
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constants). However, the off-diagonal equation E l r = implies that e = 0. Thus, we have proved the 
Birkhoff's theorem for planar transverse section, i.e. for 7 = 0, as well. It is quite remarkable that, 
even though the field equations @, are very complex in general, it admits a Birkhoff's theorem! 

4 Entropy function 

The C-function was first introduced by [13] . in the context of QFT's in 1 + 1 dimensions, who showed 
that under the RG flow to lower energies, the C-function is a monotonically increasing function of the 
couplings of the theory. At the fixed points of the flow, the C-function reaches an extremum and equals 
the central charge of the Virasoro algebra corresponding to the infinite dimensional group of conformal 
transformation in two dimensions. Later, Sahakian [J3] gave a covariant geometric expression for the 
C-function for theories which admit an holographic description. 

In [TS], Goldstein et al, gave a simple expression for the C-function for static, asymptotically flat 
solutions of Einstein's gravity in four dimensions. They showed that when coupled to matter fields 
satisfying null energy condition, this function is a monotonically increasing function of the radial 
coordinate and coincides with the entropy when evaluated at the horizon. This work was generalized 
in the context of AdS/CFT in [IB]. Recently, it was shown in [17], that C-functions also exists for 
static spherically symmetric, asymptotically flat spacetimes in Lovelock gravity. Moreover, the authors 
showed that there is a non-uniqueness in the C-function for second or higher order Lovelock theory. 
They have further shown the existence of two possible C-functions, provided the matter field satisfies 
appropriate energy conditions. Here, following the same lines of argument as in [IT], we show that one 
such C-function also exists for the theory defined^ by (pQ) . This is evident since the field equations for a 
spherically symmetric spacetime in our theory have the same functional form as that of a generic pure 
Lovelock theory [T7J. However, in five dimensions, the cubic Lovelock Lagrangian, being identically 
vanishing, does not give any field equations. Nevertheless, the theory defined by ([1]) mimics the cubic 
Lovelock theory for spherically symmetric spacetimes. 

Consider a metric of the form 

ds 2 = _ a 2 (r) dt 2 + + fe2 ^ 

a z (r) 

where S3 is a Euclidean space of constant curvature 7 = ±1,0. For 7 = 1, the spacetime (|25p has a 
spherical symmetry. The relevant components of the field equations ((4j) for this metric read 

E\ = -1 (7 - a 2 b' 2 ) 2 (7 + 3a 2 bb" + 3abb'a' - a 2 b' 2 ) (26) 

E\ = - i (7 - a 2 b' 2 ) 2 (7 + 3abb'a! - a 2 b' 2 ) (27) 

3 In a very recent paper |18] , the author constructed a cubic generalization of the BHT new massive gravity in three 
dimensions, by demanding the existence of a C-function. 
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Suppose the theory is coupled to a matter field, which satisfies the null-energy condition 



T ab ee > o (28) 

for all null- vectors £ a . This implies the following inequality: 

E\ - E\ = -^-(7 - a 2 b' 2 ) 2 b" = T\ - T\ > 0. (29) 

Now, if the metric (|25p describes a black hole, then due to cosmic censorship, b(r) / on or outside 
the horizon r = r + . Without loss of generality, we can the assume b(r) > on the horizon. For 
an asymptotically flat black hole, as r — > 00, b(r) — > ±r. First consider the case b{r) — > —r as 
r —7- 00. Since b{r) is assumed to be positive on the horizon, b{ro) = for some r+ < ro < +00, 
which is discarded by cosmic censorship. Hence, b(r) — > r as r — > 00. Now, if b(r) is not a monotonic 
function of r, then there must exist at least one local minima, i.e., b' (r c ) = with b" (r c ) < 0, for 
r + < r c < +00. However, this is ruled out by (f29l) since a 2 is positive outside the horizon. Thus the 
monotonicity of b is proved for the theory coupled to matter fields satisfying the null-energy condition. 

We now compute the entropy of a static black hole of the form ([2"o]) . using Wald's formula |T9"].j20j, 
which is given by 



f dC 

S = -2ttk / — e ab e cd i , (30) 

where e a b is the binormal to the horizon cross-section and e is the volume form induced on the spatial 
cross section S3 of the horizon at r = r + . Using the Lagrangian ([2]), we compute the curvature 
components and we obtain 



S = 127TK 



(7 - a 2 b 



2U2\2 



b 

'2 



Vol(E 3 ) , (31) 



r=r_|- 



= 12ttk— J — Vol (S 3 ) (v a = at r = r+). (32) 
b(r + ) 

Since the C-function is a function of the radial coordinate r, which matches the entropy of the black 
hole when evaluated on the horizon r = r+, one can extend the entropy formula for arbitrary r. The 
C-function is then given by 

(ry _ a(r) 2 b'(r) 2 ) 2 

C(r) = 12itk U { ' I { 1 ' Vol (S 3 ) (33) 



Let us now check the monotonicity of f|33|) as a function of outward radial coordinate, following along 
the lines of Ref. [17J. First note that using the field equation (126p . one can write tt component of the 
stress energy tensor as 

Tu = P U\X - ^)1 , (34) 
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where X := 7 — a 2 b' 2 . Now, differentiating C{r) we obtain, 

1 27T 

C'ir) = -^K(2bXX' - X 2 b')Vol (S 3 ) 



^(^-^(X-^))K 0!(S 3) 

) UoZ (S 3 ) . (35) 



46'vr / 2 4& 6 T t ^ 

-K -A 



b 2 \ kXo 2 

Using null energy condition, we had shown that b' > (See the paragraph below Eq. (j29[) ) and since 
weak energy condition implies Tu > 0, then for both 7 = 0,-1, X < 0, and hence from Eq. (|35p . 
we see that C > 0. This proves that, for 7 = 0, —1, the function C(r) is a monotonically increasing 
function. However, our analysis is inconclusive for 7=1. 



5 Black hole thermodynamics 

In this section, we explore the thermodynamics of the black hole (|13p . which as proved in the previous 
sections, is the unique solution for the theory ([1]) within the family of spacetimes (|14|) . 
The temperature of the black hole (|13p is given by 

, / / \ 2/3 x ' 
1 / / r x ' 



47T \ V r A ' 



r=r+ 



, (36) 



which is also the case for spherically symmetric black holes in pure Lovelock theories for arbitrary 

2 



order k < ^J- . Using Wald's entropy ([32]) one obtains 



1 2-7T 

S = nVol (S 3 ) . (37) 

r+ 

And asuming the validity of the first law dM = TdS, one finds that the mass of the black holes is 
given bjfl 

M = V ° l P — , (38) 

where we have fixed the integration constant Mq in such a way that when r + — > 00 (i.e. flat space) 
the mass vanishes. 

Since the mass in terms of the temperature is given by 

M = 36tt 2 T 2 kVoI (S 3 ) , (39) 



4 It would be interesting to compute the mass by more standard methods like the Hamiltonian analysis. 
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the specific heat C = dM/dT is positive 



C = 72tt 2 kT Vol (£ 3 ) , (40) 
which implies that the black hole is thermodynamically locally stable. 

6 Nonhomogeneous combinations: Asymptotically AdS black holes 

In five dimensions the most general Lagrangian giving rise to second order field equations is given by 
an arbitrary linear combination of the Gauss-Bonnet, the Einstein's and cosmological terms. We now 
look for nontrivial, spherically symmetric solutions when the cubic Lagrangian ^) is supplemented 
by a linear combination of the above terms. We find a new asymptotically AdS black hole for a 
particular combination, which in addition to the event horizon has a Cauchy horizon. The field 
equations obtained for the theory considered, are given by 

£^ u := c 3 E^ u + c 2 GB^ v + ciG^ u + cog^ u = , (41) 

where is defined in Eq. @, G^ v is the Einstein's tensor and the Gauss-Bonnet term is defined by 

GB^ := 2RR 11V - AR pp R p „ - 4R S p R p ^ + 2R^R V & - ^(RpS^R^ ~ 4R pS R p5 + R 2 ) • (42) 

Equation (I4ip describes the most general cubic theory in five dimensions [6], whose field equations are 
of second order, for static spherically symmetric spacetimes ([6]). 
The constant curvature solutions of this theory 

Kb = A (m - W) , (43) 



J a u B w a w /3 

fulfill 

2c 3 A 3 - 12c 2 A 2 - 6ciA + c = . (44) 

Generically, there are three constant curvature solutions with different radii (inverse of different cos- 
mological constants), describing three different maximally symmetric spacetimes, which corresponds 
to (A)dS or flat space, depending on whether A is (negative)positive or zero, respectively. In analogy 
with Lovelock theories, it is natural to expect that, the space of the solutions is enlarged when the 
three different vacua of the theory degenerate into one [11], |12] , [21j-|28j, which occurs when 



c 2 = — - , and c 3 = -4-i . (45) 
co eg 



In such a case Eq. (|44|) factorizes as 

(c -2c 1 \f 



. (46) 
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Consequently for cq 7^ 0, one obtains 



A = S- (47> 



Assuming for simplicity / (r) =5 (r) in ([6]) we integrated the field equations to obtain the following 
solution 

ds 2 =-(£- cr 2 " + 7) di 2 + r3 ^ + r 2 ^ (4g) 

V / 73- - cr 2 / 3 + 7 

where c is an integration constant, 7 is the curvature of S3 and I 2 := — — is the squared AdS 
radius which is assumed to be positive. For c = 0, the spacetime is locally AdS, and in this case for 
7 = — 1, the metric reduces to the massless topological black hole [29]. For 7 = 1, the metric (|48p 
is asymptotically AdS with a slower fall-off as compared with the Henneaux-Teitelboim asymptotic 
behavior [30]. This again, is similar to what occurs in Lovelock theories [31j,[32j. The spacetime (|48p 
has a curvature singularity at the origin which could be covered by one or two horizons depending on 
the values of c and 7. 

For 7 = 1, and —00 < c < 6 (p) 1 ^ 3 , the metric ([4*5]) describes a naked singularity. In the case 

1 /3 

c = 6 (p-) the spacetime ()48|) describes an extremal black hole with a degenerate horizon located 
at r + = r_ = 21. In the range c > 6 (p-)^ 3 the metric (j4~8|) has an event and a Cauchy horizon, which 
cover the timelike singularity at the origin. 

For vanishing 7, the singularity at the origin becomes null, and for positive c this singularity is 
hidden by an event horizon located at r + = (c/ 2 ) 3//4 . 

Finally, in the case 7 = — 1, there exists an event horizon at r + for any value of c, which covers a 
spacelike singularity at the origin. 



7 Generalization to arbitrary higher order 

Some aspects of the theory defined here, are common to other, well behaved, theories of gravity, 
such as Lovelock theory and BHT new massive gravitjjf] where the trace of the field equations is of 
second order. Furthermore, there is an interesting similarity with the cubic Lovelock theory, which 
was exploited in the construction of the C-function. The "spherically symmetric" solution of the pure 
cubic Lovelock theory, is given by [38] 

l + — )dt 2 + - + r 2 6Z\ D _ 2 , (49) 

where c is an integration constant. This is valid for D > 6. Nevertheless, if one insists on considering 
D = 5 in (|49p . one obtains the metric (jlip . This is also the case for pure BHT new massive gravity 
[39j , [41] , when the spherically symmetric solution to pure Gauss-Bonnet field equations is "extended" 
to D = 3 [31]. The thermodynamics of the black holes found here, also reveals some similarities with 
the ones of pure BHT, where the specific heat is positive (and linear in the temperature), implying 

5 Note that for the pure BHT theory there is no Birkhoff 's theorem which is explicit from the existence of gravitational 
solitons [39| . This non-uniqueness further allows the existence of a very interesting Lifshitz black hole [40] . 
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the thermal stability of these black holes. Naturally, due to the remarkable resemblence between the 
BHT new massive gravity in three dimensions and our theory in five dimensions, one can't help but 
wonder if there exists suitable generalizations to arbitrary higher order. We show, in this section, that 
the answer is affirmative by presenting a recipe to construct generalizations of our theory to arbitrary 
higher order. In Appendix B, we give explicit results for the quartic case. 

First, let us recall how the quadratic invariant K := 4R ab R a t, — n^z^ R 2 can be constructed, which 
in D = 3 serves as the Lagrangian for the BHT new massive gravity. The key step here is to realize 
the following identity in arbitrary dimensions. 

^abcdri c i I ^ ~ ^\ I a r>ab r- ^ 



C abcd C abcd = S 4 +[ — {4R ab Rab - Jd^I) R 1 ' m 



which can be rewritten as 



D-l \D-3 



Cabcd s~i c 



(51) 



where £4 := R 2 — 4R a bR ab + R a b c dR abcd is the four dimensional Euler density and C a bcd is the Weyl 
tensor. At first sight, it seems that in three dimensions, the right hand side takes a 0/0 form since 
both the Weyl tensor and the four-dimensional Euler density vanishes identically in D = 3. However, 
if one expands the right hand side in terms of the Riemann curvature tensor then it factorizes by 
{D — 3) which cancels the one in the denominator of the preceeding factor and we are left with the 
combination on the left hand side. Note that the difference of the quadratic Weyl invariant and the 
four-dimensional Euler density can also be written as 

}_^a 1 bia 2 b 2 ( (~<c\d\ (~ic 2 d 2 _ ryc\d\ -nC^dzS (rn\ 
22°c 1 d 1 c 2 d 2 ^0161^0262 ^0161 -"-0262 J { > 

where 5..'. is the totally atisymmetric tensor. 

Now, we are ready to generalize the identity (|5ip for higher oder. First, consider the following 
invariant of order k 

J_ eai&i---a fe 6 fc ( f~ic\d\ /~<c k d k _ pcidi nc k d k \ /t;o\ 

2 fc °cidi-c k d k ^ U aibi ' ' ' ^a k b k ^aibi " " " n a k b k ) l ^ 

Obviously, the above invariant vanishes in dimensions lower than 2k. However, if one expands the 
Weyl tensor in terms of the Riemann tensor, then it can be factorized by (D — 2k + 1). This can be 
seen as follows. Consider the basis set of A:-th order Riemann invariants in arbitrary dimensions. In 
D = 2k — 1, not all elements of this set are linearly independent. In fact, the basis set contains one less 
invariant than in D > 2k. This is beacuse of the vanishing of the fc-th order Lovelock density. Now, 
after the expanding in terms of the Riemann tensors, the term (|53|) will not contain any (Riemann)^. 
So, this invariant cannot vanish identically in D = 2k — 1 unless it is factorized by {D — 2k + 1 



ffl 



6 This argument cannot be extended to dimensions 2k — 2 since one obtains another identity involving the Riemann 
invariants which is obtained by contracting the Ricci tensor with the (k — l)-th order Lovelock equation. 

7 Further expanding all the Weyl tensors, one can convince one self that the dimensional dependence of the coeficient 
of the term with k — 1 Riemann tensors and one Ricci tensor must be (D — 2k + 1) / (D — 2). 
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We can now divide this factor out to get a non-vanishing invariant in D = 2k — 1. Thus, we write the 
kth order generalization of K by evaluating 



1 ( 1 



f.a\b\---a k b k ( (~ic\d\ (~<c k d k _ pCidi nc k d k \ (KA\ 

2 k \ D - 2k + I J c i rf i- c fc rf fc V aifel ' ' ' a ^ k "0161 " " " U a k b k J K^) 

in D = 2k — 1. Now, by construction, the trace of the field equation arising from the above invariant is 
of second order in all dimensions. Moreover, for static spherically symmetric spacetimes this invariant 
must be a sum of the kth Lovelock invariant and a term proportional to (Weyl) k . This is because for 
static spherically symmetric spacetimes, all the Weyl invariants are proportional to each other |43j . 
Thus, if one subtracts an appropriate multiple of (Weyl) k from this invariant, then all the components 
of the field equation, arising from the resulting invariant will be of second order. Since, for k > 4, 
there are more than one linearly independent Weyl invariants in dimension 2k — 1, one can do this in 
several ways. One convenient choice is 

r . ^ 1 I xaibi—akbk ( (~iC\d\ r~ic k d k _ pcidi jjc k d k \ _ /~ia k b k s~ia\b\ ^a k _ 

2 k \D - 2k + l) c idi - c fc d fc \^ai&i ' ' ' ^a k b k ' ' ' n a k b h J "^aidi ^a 2 b 2 ' ' ' ^a k b. 



(55) 



where 



_ (D - 4)! [k(k - 2)D(D - 3) + k(k + !)(£> - 3) + (D - 2k)(D -2k- 1)] 
ak ~ (D-2k + l)\ [(D-3) k ~ 1 (D-2) k - 1 +2 k - 1 -2(3-D) k - 1 } { } 

Again, note that the above invariant vanishes identically in D < 2k — 2, whereas in dimensions 
D > 2k, it can be expressed as a linear combination of the Weyl invariants and the 2fc-dimensional 
Euler density. After replacing the ansatz 

ds l = _at (r) / (r) dt 2 + — - + r 2 cffi2 2fc _ 3 (57) 

in (|55p and carrying out the variation of the action with respect to / (r) and iV (r), one respectively 
obtains 

( 7 -/) fc - 1 JV' = 0, (58) 
(7 " ft' 1 [(D - 2k - 1) (7 - /) - krf] = . (59) 

The non-trivial branch of solutions for D = 2k — 1, is given by 

ds 2 = - (cr^ + 7) dt 3 + + r 2 dSL_3 , (60) 

where c is a integration constant and dE( 2 fc_3) is the line element of a Euclidean (2k — 3)-dimensional 
space of constant curvature % 7 = ±1, 0. For positive c and 7 = —1, this describes a topological black 
hole with a horizon located at r = r + = c~? . The temperature of the black hole is l/27rkr + . One can 



In three dimensions for BHT, 7 is an integration constant. 
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compute the entropy using Wald's formula and then obtain the mass with respect to the locally flat 
background by applying the first law, 



S oc 



Vol (S 2 fc-3) 



M oc 




(61) 



8 Further comments 

In this paper we have investigated a new interesting theory of gravity, which is cubic in curvature. 
This theory as we have shown, has several remarkable characteristics such as having second order field 
equations for static spherically symmetric ansatz, admittance of Birkhoff 's theorem and existence of 
a C-function for the black hole solution. This theory is the unique cubic theory in five dimensions, 
for which the field equations for the static spherically symmetric spacetimes are of second order [6]. 
As in the case of Lovelock theory, the admittance of Birkhoff 's theorem [9], |10| . suggests the lack of 
the spin-0 mode in the linearized theory [33], [53]. A definite confirmation to this assertion requires 
a full Hamiltonian analysis, which can be performed for example along the lines of Ref. [35J , or in a 
spherically symmetric minisuperspace approach [36], which is straightforward due to the second order 
nature of the theory in this setup. This can be seen from the Lagrangian, where all the terms that 
are second order derivatives in the metric functions are of the form H{q) q, which can be integrated 
by parts to obtain a first order Lagrangian. 

It is natural to expect that, when perturbed around flat space (up to the leading order), our theory 
will possess ghost degrees of freedom, since generic perturbations will break the spherical symmetry 
and involve fourth order derivatives. Nevertheless, since this assertion is background dependent, it 
would be nice to look for a ghost-free background, in analogy with Topologically Massive Gravity [57] . 

In the nonhomogenous combination, due to the presence of additional scales, it is natural to expect 
that the black hole (|48p will have different phases depending on the sign of the specific heat, as is the 
case for the black holes studied in [31]. It will be also interesting to prove Birkhoff's theorem for all 
the higher order generalizations. Work along these lines is in progress. 

It would be interesting to explore the dimensional reduction of the nonhomogeneous theory to four 
dimensions, along the lines of Ref. [42j . 

Note: After the first version of this work was submitted to arXiv, another paper [33J appeared 
where the same cubic invariant in five dimensions is presented. The authors generalize the cubic 
invariant in higher dimensions by requiring second order field equations for static spherically symmetric 
spacetimes. However, their generalization in higher dimensions is nothing but a particular linear 
combination of the six-dimensional Euler density and the two linearly independent Weyl invariants. 
As mentioned above, since for static spherically symmetric spacetimes, both the Weyl invariants are 
proportional to each other, addition of the two invariants with a particular choice of the relative factor 
does not contribute to the field equations. In this case, one obviously obtains the field equations for 
the cubic Lovelock theory. 
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A General static, spherically symmetric solution for the 
non-homogeneous combination 

In this appendix, we consider the general static, spherically symmetric solution for a generic non- 
homogeneous combination (|41|) . For simplicity, we restrict to spacetime metrics of the form @ , with 
/ = g. In this case the components of the field equations £ t and £ r r are equal and reduce to 

£\ = £\ := 2c r 3 - 3 C1 (r 2 ( 7 -/))'- 6c 2 (( 7 - /) 2 )' + c 3 (r~ 2 ( 7 - /) 3 )' = , 
which, after defining F := 7 — / reduces to 

2c r 3 - 3ci (r 2 F)' - 6c 2 (F 2 )' + c 3 (r~ 2 F 3 )' = . 

This equation can be trivially integrated to obtain the following algebraic equation 

c r 6 - 6cir 4 F - 12r 2 c 2 F 2 + c 3 F 3 + 2{ir 2 = , (62) 

\x being the integration constant. The field equations with indices along the manifold S3, i.e. • are 
given by 

*J = ^2 (^r)' «i ■ (63) 

Therefore, the solutions to equation (|62p . trivially solve equation (|63p . This again is analogous to Love- 
lock theories, in which the problem reduces to solving an algebraic equation, given by the Wheeler's 
polynomial [15]. Generically there are three branches, which are asymptotically locally a spacetime of 
constant curvature Ai, A 2 or A3, when the equation (|62p has real roots. 

Now, let us examine a particular case for which the solution takes a simple form. This is the 
case when the pure cubic theory is appended by a cosmological term, i.e. c\ = and c 2 = in (|4"Tj) . 
This is the simplest case which admits an asymptotically locally AdS solution. In this case, the cubic 
equation reduces to 

c r 6 + c 3 F 3 + 2[ir 2 = , (64) 
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which is solved by 



/W=7+^(l-^) 1/3 , (65) 



where we have defined the AdS radius as I 2 := (c^/cq) 1 ^ 3 (chosen to be positive) and the integration 
constant \x has been replaced by fl = — / (3cs). The spacetime described by ([65]) is asymptotically 
locally AdS, and describes a topological black hole for a certain range of the parameters. Note 
that expanding around infinity, the subleading term goes as a Schwarzschild-Tangherlini term (/2/r 2 ), 
suggesting \x as the mass parameter. 

B A new quartic theory of gravity in seven dimensions 

Here, we present the generalization of our theory to quartic Lagrangians in seven dimensions. Consider 
the following basis of quartic invariants 



T — TfPqbs p a u n v w p r _ npqbs n a u n v w n j _ npqbs n au n V w n 

— - n ' - n 'p b q JrL uvswi ±J 2 — J J^p b -<</ a u ^qvswi ^3 — J - n 'pg a -^svuwi 

T — npqbs n au n vw p r _ npqbs p au p vw p r _ npqbs p a puuiw p 

J-14. — JX ^pq J^ba ^suvwi J -'5 — I '- ^pq ^au ^bsvwi ^6 — ^pqb s^'avwai 

T — ( npqbs p \ 2 r _ ppg nbsau p u p r _ ppg nbsau p u p 

A/7 — I it n-pqbs I j -kg — -ft -ft -ft;, ap rL svuq , -Lg — JX JX JXjj S pJXauvqi 

T — DPI b s nauv p r _ p npqbs p a u p r _ p npqbs p au p 

Mo — -K -ftp g b-Kauvs, Ml — -K-K ftp b Mqasu, J->\2 — tVx tipq Mbsau, 

t ppg pfes pa « p r ppg p6s pa « p r ppg pfes paw p 

-^13 — -ft K H p b-Kaqus, -^14 — Jx Jx JX p q Jx a bus, -^15 — -ft -ft -ft pb^auqs, 

t ppg p b nsau p r ppg p nbsau p r p ppg pfcsa p 

M6 — -ftp -K q-tisaub, ^17 — -K -ftpg-ft -Kfosan, -^18 — -K-K -ft p rlbsaq, 

L19 = R 2 R pqbs Rpqb s , L20 = R pq R bs Rb a R P sqai L21 = RR pq R bs Rpb qs , 

L22 = R pq R p b R q s Rbs , L23 = (R pq R P q) 2 , L24 = RR pq R p b Rqb, £25 = R 2 R pq R P q, £26 = -ft 4 - (66) 

In seven dimensions, they are not linearly independent, since they are related to the eight-dimensional 
Euler density 

-B 8 := -96Li + 48L 2 - 96L 3 + 24L 4 + I8L5 - 48L 6 + 3L 7 + 384L 8 - 192L 9 + 192L 10 - 32L n + I6L12 
+ 192Li3 - 192L M + 96Li 5 + 192Li 6 - 24L i7 - 96Li 8 + 6L 19 - 384L 20 + 96L 2 i - 96L 22 + 48L 23 
+ 64L 24 - 24L 25 + L 26 . 

which vanishes identically in dimensions lower than eight. Now, we consider a Lagrangian, constructed 
by taking a linear combination of all the invariants from the set (|66p . with the following choice of 
coefficients 

:=16Li + 38L 2 + - 14L 6 - ^L 7 + 16L 8 - 32L 9 + 2L 12 - 24L 14 + 8L15 + 16L 16 + ^-L 17 

2 lb 1U 

61 8 153 16 121 57 

" ^^19 + 5^22 " ^^23 " y^24 + ^^25 " M ^26- (67) 
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In analogy with the quadratic (BHT new massive gravity) and the cubic case, the trace of the field 
equations obtained from the above Lagrangian, being proportional to the Lagrangian itself, is of second 
order. In fact, there are seven Weyl invariants, constructed by taking complete contractions of four 
Weyl tensors, which also has this property in arbitrary dimensions. In terms of the basis of invariants 
([66]) . they are given as 

8/ L 9 \ 4L 10 2 /_ L 12 \ 22L 13 , 16L U 4L 15 12L 16 2L 17 
Wr-Lx - - [L 8 - — — + — (L u — — - + — — — — + 



5 V 4 J 5 15 V 4 J 25 25 25 25 25 

8Lis Li 9 172L 20 _ 42L 2i 201L 22 _ 103L 23 _ 416L 24 52L 25 _ 13L 26 
25 50 125 125 + 625 625 1875 + 625 3750 

w T , 8/ L 9 \ 2 / L 12 \ 12L 13 28L 14 4L 15 8L 16 4L 18 L 19 

Wo:=L 2 H — ^8 -^11 H 

5 V 4 y 15 V 4 / 25 25 25 25 25 150 

136L 20 46L 2 i 138L 22 104L 23 308L 24 148L 25 37L 26 

125 + 125 625 + 625 + 1875 1875 11250 

8L 8 2L 9 2Lio 2Ln Li 2 12Li 4 12 / L15 \ 3L15 2Li6 L17 

W3:=L 3 L13 

5 5 5 15 30 25 25 V 2 7 5 5 25 

6L 18 L 19 144L 20 34L 21 152L 22 76L 23 3 32L 24 122L 25 61L 26 

~25 75~ 125 125~ 625 625 1875 1875 ~ 22500 

_ 8L 9 2L 12 16/ Li 5 \ 24L 15 , 16L 16 8L 18 L 19 96L 20 , 16L 21 

Wa:=La Li 3 H 

5 15 25 V 2 J 25 25 25 75 125 125 

168L 22 24L 23 96L 24 68L 25 17L 26 

625 625 625 1875 11250 

_ 16L 9 4L 12 32L 13 32L 15 32L 16 16L 18 2L 19 192L 20 , 32L 21 336L 22 

W^:=Lk 

b b 5 15 25 25 25 25 75 125 125 625 

48L 23 192L 24 136L 25 17L 26 

625 625 1875 5625 
w . T 8L 10 16L 14 28L 16 4L 17 8L 18 2L 19 112L 20 32L 21 196L 22 108L 23 

Wa'.=LF. 

5 25 25 25 25 75 125 125 625 625 

112L 24 136L 25 17L 26 



TTT T 8-L17 2-Lig 16-L 2 ':i 8-L95 -^2fi / n 

W 7 :=L 7 ii + — + fl _ + _f°. (68) 



625 1875 5625 

15 25 75 ' 225' 
Note that the Weyl invariants satisfy the following identity 

E 8 := -96W X + 48W 2 - 96W 3 + 24iy 4 + I8W5 - 48W 6 + 3W 7 = 0. (69) 

This implies that only six of the above Wi's are linearly independent in seven dimensions. However, 
the invariant (|67p is linearly independent of the set of Weyl invariants fj68[) and consequently does not 
transform covariantly under Weyl rescalings. 
Let us consider the following action 

I 4 = k 4 / ^gL^d 7 x. (70) 
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Now, consider a static spherically symmetric spacetime described by the line element 

ds 2 = -N(r)f(r)dt 2 + _ - + r 2 d£ 2 (71) 



where CE5 ~ is the line element of an Euclidean five-dimensional space of constant curvature 7 = ±1, 0. 
Using the minisuperspace method [36], we obtain the following two second order equations for f(r) 
and N(r): 

N' 3 = 0, (72) 
( 7 - f(r)f[2rf'(r) + 7 - f(r)]) = 0. (73) 

The nontrivial solution for the pure quartic case is given by: 

lb .__( 1T v. +7 )tf + _*L_ + ^, 

where c is an integration constant. As in the pure cubic theory, for the topological case 7 = — 1 
and positive c, this metric describes an asymptotically locally flat black hole with horizon radius 
r = r + := 1/c 2 , that can be rewritten as 

ds 2 = -((-) - l) dt 2 + - + r 2 d^ 2 . (74) 



In this case, the horizon is described by a quotient of the five-dimensional hyperbolic space H5/T, 
where T is a freely acting discrete subgroup of O (5, 1). 
The temperature of the black hole is 

T- 1 1 

8-7T r + 

Using Wald's formula for the entropy and then the first law to compute the mass, we obtain 

S = ^9l K4 Vol (£ 5 ) (75) 

respectively. Finally, the specific heat of this black hole is given by 

C = 3840vr 2 K 4 r Vol (S 5 ) , (77) 

which, being positive, implies that the black hole is thermally stable. Note that the functional depen- 
dence of all the expressions remain the same as their five-dimensional, cubic counterpart. 

Let us further consider a generic linear combination of the quartic and all possible lower order 
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Lovelock terms: 



where L° := 1, L 1 := R and 



r 4 

I 7 = k 4 y/^g^OiL^tfx, (78) 

i=0 



L 2 :=R abcd R abcd - 4R ab Rab + R 2 , (79) 
L 3 :=2R abcd R cdef R ef ab + m ah cd R ce bf R d f ae + 2AR abcd R cdbe R e a + (80) 

o p jyabcd p _i_ OA oabcd p p i ic P"?> P pc ion pafe p . p3 

oKK tx abcd + Z4/t H ac K bd + 10/1 ilbc-fl a — 1ZKH H ab + /I . (51 J 

The field equations for spherically symmetric ansatz (|7ip reduce to 

- r 5 a + 5 (r 4 F)' ai - 60 (r 2 F 2 )' a 2 + 120 (F 3 )' a 3 - 30 (r" 2 F 4 )' a 4 = , (82) 

where F = j — f . The above equation can be trivially integrated as 

-r 8 a + 30r 6 F«i - 360r 4 F 2 a 2 + 720F 3 r 2 a 3 - 180F 4 a 4 + 6r 2 /i = . 

Here fi is an integration constant. Obviously, depending on the roots of the above polynomial equation, 
the metric may then describe a black hole. 

The maximally symmetric (see f)43|) ) solutions of this theory fulfill 

a + 30Aqi + 360A 2 a 2 + 720A 3 a 3 + 180A 4 a 4 = . (83) 

Generically, a spherically symmetric solution will asymptotically approach a maximally symmetric 
background of constant curvature Ai, where Aj is a real root of (|83l) . When the coupling constants a% 
are such that 

4 3 2 

ce = 180Hr, a x = 24^|, and a 2 = 3^ , (84) 
ai ai a 4 



equation (f8"3"j) reduces to 

(a 4 + a 3 A) 4 



, (85) 
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d* 2 = - ( F - ^r 1/2 + 7 J ^ 2 + r2 _ z 1/2 , : + ^ 2 ^i, 7 , (86) 



and the four maximally symmetric vaccua of the theory coincide. The spherically symmetric solution 
is then given by 

V 2 _ 1/Q \ ,. 9 dr 2 

^ - fir 1 /^ + 7 

where I 2 := a^/a^ is the square of the AdS radius and jl is the rescaled integration constant. Depending 
on the values of the parameters, the metric may then describe an asymptotically AdS black hole or 
a topological black hole. The fall-off at infinity is slower than the one in GR and the spectrum of 
spacetimes is the same as in the cubic case, which we described at the end of section ([6j). However, as 
mentioned earlier, for static spherically symmetric spacetimes, all the Weyl invariants are proportional 
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to each other. Specifically, in seven dimensions they are related by 



w 1= w 1 = w 1= w ±= w^ = w 1= w 1 _ 

93 191 11 226 452 328 1000 ' 1 ' 

This means that one can always add an arbitrary combination of the Weyl invariants YlJ=i C"iW{ to 
the Lagrangian ([67|) . such that 

93Ci + 191C 2 + HC 3 + 226C 4 + 452C 5 + 328C 6 + IOOOC7 = 0. (88) 

without affecting the field equations and their spherically symmetric solutions. 
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